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Abstract 

The exact solution of directed self-avoiding walks confined to a slit of finite width and interacting 
with the walls of the slit via an attractive potential has been calculated recently. The walks can be 
considered to model the polymer-induced steric stabilisation and sensitised floculation of colloidal 
dispersions. The large width asymptotics led to a phase diagram different to that of a polymer 
attached to, and attracted to, a single wall. The question that arises is: can one interpolate between 
the single wall and two wall cases? 

In this paper we calculate the exact scaling functions for the partition function by considering 
the two variable asymptotics of the partition function for simultaneous large length and large 
width. Consequently, we find the scaling functions for the force induced by the polymer on the 
walls. We find that these scaling functions are given by elliptic ^-functions. In some parts of the 
phase diagram there is more a complex crossover between the single wall and two wall cases and 
we elucidate how this happens. 

1 Introduction 

The problem of a single polymer confined between two walls and interacting with those walls has been 
considered for at least 35 years pp. One reason for this is its use as a model of the stabilization of 
colloidal dispersions by adsorbed polymers (steric stabilization) and the destabilization when the poly- 
mer can adsorb on surfaces of different colloidal particles (sensitized fiocculation) . Until recently, even 
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when one substitutes directed walks for the more canonical self-avoiding walks in a two-dimensional 
lattice model of this phenomenon, the only cases to be considered exactly have been special cases 
where the interaction with the two surfaces are equal. Recently, Brak et al. [2] have calculated the 
generating functions for a directed self-avoiding walk confined by two horizontal walls on the square 
lattice where separate Boltzmann weights a and b were associated with visits to the lower wall and 
upper wall respectively, with various restrictions on the end-points of the walks. 

The dominant singularity of the generating function of any of the subcases considered by Brak et 
al. [2] leads to the calculation of the free energy and the force induced by the polymer to the walls. By 
considering the dominant singularity one is effectively considering the infinite walk length limit. The 
dominant singularities of the generating functions were analysed asymptotically for large widths. In 
the infinite width limit a novel phase diagram was obtained, different to the one obtained by analysing 
a polymer in a half-plane geometry. In Figure [T] the phase diagram obtained by Brak et al. [2] is 
given. For small a and b the polymer is desorbed from both walls while for large a, respectively large 
b, the polymer is absorbed onto one wall or the other. It is interesting to note that new numerical 
results 0] and rigorous results [5] show that undirected self-avoiding walks demonstrate very similar 
behaviour to the exact solution of the directed walk model. In the infinite width limit, referred to as 
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Figure 1: Phase diagram of the infinite strip where the separation of the walls is made large after the 
limit of infinite walk length is taken. There are 3 phases: desorbed (des), adsorbed onto the bottom wall 
(ads bottom) and adsorbed onto the top (ads top). Our notation for the various regions and transition 
lines are marked. The transition lines at a = 2,6 < 2 and b = 2, a < 2 marking the boundary of the 
desorbed region are second order phase transitions with a jump in the specific heat on crossing the 
line while the line marking the boundary at a = b,a > 2 of the two adsorbed regions is a first order 
transition. 
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the infinite slit (in two dimensions), the (reduced) free energy «(a, b) is given by 



K is (a,b) 
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In the half-plane geometry the limit of infinite wall separation is effectively taken before the limit of 
infinite walk length (see Figure Ej). The results in [2j imply, unusually, that the order of the two limits, 
polymer length to infinity and wall separation to infinity, are not interchangable. In fact, the free 
energy depends only upon the value of a and is given by 
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Hence, for b > 2 and a < b (denoted region T in Figured]) the infinite slit and half-plane free energies 
are different. The question that naturally arises is whether there is a scaling function that interpolates 
between these two limits and whether this extends to the region where the free energies differ. To do 
this one needs to consider the finite length partition functions rather than the generating functions. 
An exact expression for the finite length partition function has now been calculated [6|. However, it is 
not easy to see how one could analyse this expression asymptotically, especially for large separations. 
Returning to the undirected model, in three dimensions, a scaling theory [3] valid for the desorbed 
region T> of the infinite slit and its boundaries has been shown (numerically) to hold. 



des 



ads 



a. 



a 



Figure 2: Phase diagram of the half-plane (one wall) problem where the separation of the walls is made 
large before the limit of infinite walk length is taken. There are 2 phases: desorbed (des) and adsorbed 
onto the bottom wall (ads bottom). The boundary of the two phase is a second order phase transition. 



In this paper we calculate the two variable asymptotics for large wall separation and large polymer 
length of the partition function of one of the directed walk models considered by Brak et al. [2] . The 
paper is set out as follows. In Section 2 we set the stage for our calculations: we recall the exact 
definition of the model and the calculated generating function, and introduce the notation which we 
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will use in the remainder of the paper, in particular, a parametrisation in more convenient variables. 
Much of the analysis hinges on the understanding of the singularities of the generating function; this 
will be discussed in Section 3. In Section 4 we show how the contour integral expression for the partition 
function can be reformulated in terms of residues of the singularities of the generating function. We 
discuss special cases for specific values of a and b, where one obtains simple expressions for the partition 
function, and also give an exact expression for general a and b. In Section 5 we calculate the scaling 
function of the partition function for the special cases, and then derive the scaling function in the 
various regions of the phase diagram for general a and b. In Section 6 we discuss the results in the 
light of finite size scaling theory and, in particular, discuss how in each region of the phase diagram the 
scaling results interpolate between the single and double wall models. Importantly, we demonstrate 
that in the desorbed region of the infinite slit and on its boundaries the scaling theory proposed for 
undirected SAW in a three-dimensional slab holds exactly for our directed model, with appropropiate 
exponent substitutions. 



2 The model 

Brak et al. [2] considered three different end-point restrictions. In this paper we shall restrict ourselves 
to the case when both end-points are attached to the same wall: in [2] these were referred to as loops 
(see Figure [3]) . If £™ is the set of loops of fixed length n edges in the slit of width w then the partition 




Figure 3: An example of a directed path which is a loop: both ends of the walk are fixed to be on the 
bottom wall. A Boltzmann weight a is associated with visits to the bottom wall (excluding the first) 
and a Boltzmann weight b is associated with visits to the top wall. 



function of loops is defined as 

Z n>w (a,b)= Y, fftovV* (2.1) 

where u(p) and v(p) are the number of vertices in the line y = (excluding the zeroth vertex) and the 

number of vertices in the line y = w, respectively. The generating function L w {z, a, b) is then given by 

oo 

L w (z,a,b) = J2z ntW (a,b) z n . (2.2) 

n=0 
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The force T ntW {a, b) is defined as 



T n w (a, b) = ^(Z n , w+ i(a,b))-log(Z n , w (a,b)) 

n 

but can be estimated from an asymptotic expression for large w as 

T - \ dZn > w (2 4) 

•r n,w — ry a \^-^J 

The aim of this paper is to extract the asymptotics of the finite size partition function Z n>w (a,b) 
by inverting equation f|2.2[) . Hence we have 

1 f dz 
Z n ,w(a, b) = — j> L w (z, a, b)—^ , (2.5) 

where the generating function L w (z, a, b) has been calculated in [2] and is 

r ( h s_ (l + q)[(l + q-bq)-(l + q-b)q w ] 

W{ ' ' ' ~ (1 + q - aq)(l + q - bq) - (1 + q - o)(l + q - b)q™ ' ^ 

with z = ■\/q/(l + q). The problem is to evaluate the above contour integral for large but finite n and 
w. 

Before we enter the calculations let us introduce some notation to keep track of the different 
regions and transitions in the phase diagram of the infinite slit (Figure [1]) . We restrict our discussion 
to a, b > 1 and label the desorbed region, a < 2, b < 2 as T>, the region where the polymer adsorbs 
onto the bottom wall, a > 2, a > b as B, and the region where the polymer adsorbs onto the top wall, 
b > 2, b > a as T. The boundary where the region T> meets region B, a = 2, b < 2, that is, when 
the polymer is critically adsorbing onto the botttom wall, we denote as Ag. The boundary where the 
region T> meets region T, b = 2, a < 2, that is, when the polymer is critically adsorbing onto the top 
wall, we denote as A3-. The boundary where the region B meets region T , a = b, a > 2, that is, when 
the polymer is equally adsorbed onto the both walls, we denote as Stb- Finally, the point where all 
three regions and all three lines meet at a = b = 2 is denoted A-rg. 

Mathematically it is advantageous to re-parametrise L w by introducing a = 1 + A 2 , 6 = 1 + fj? 
with A, \x > and q = p 2 . Hence z = p/(l + p 2 ). In what follows, we will thus work with 

L W [Z, a, 0) - L w [p, A, Li) - (1 _ AV)(1 _ _ (A2 _ p2) ^ 2 _ p2)p2w , {I. I) 

so that equation (|2.5p becomes 

Z n>w {a,b) = ±-Jc w (p,X,Lt)(l-p 2 )(l+p 2 r- 1 p % . (2.8) 

3 Singularities 

Of crucial importance for the understanding of the structure of the generating function are its singu- 
larities, i.e. the zeros of the denominator polynomial 

D w (p) = (1 - AV)(1 - MV) - (A 2 -p 2 )(li 2 -p 2 )p 2w . (3.1) 
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It is convenient to look at some special cases first, where D w {p) simplifies considerably. We have 
D W (P) = 



1 - p 4+2w for A = 0, n = , 

(1 -p 2 )(l + p 2+2w ) for A = 0, n = 1 or A = 1, p = , (3.2) 

(1 - A 2 p 2 )(l - A" V)(l - P 2w ) for Xp = 1 • 



For these special cases, all zeros are simple with the exception of the case A = 1 and p = 1, in which 
case p = ±1 is a multiple zero. For general values of A and /i, we have the following result on the 
multiplicity of zeros. 

Lemma 1. If (A,/i) / (1, 1), the polynomial 

D w (p) = (1 - AV)(1 - M V) - (A 2 -p 2 )(/x 2 -pV" (3-3) 
/ias simple zeros except possibly at p = ±1 /or a single value of w. 

Proof. The polynomial D w {p) is simply related to the orthogonal polynomials, P w (z), defined in [2]. 
In particular 

A*(p) = (i -p 2 )(i +p 2 ) w+1 p w {jT?) ■ (3 ' 4) 

Note that P w {z) is a polynomial of degree w + 1 in z, so the above expression is indeed polynomial in 
p. The zeros, pj, of D w (p) are therefore either p = ±1 or images of the zeros, z%, of P w (z) given by 



l±Jl-4zt 



A standard result on orthogonal polynomials (see Theorem 5.4.1 in [7] for example) implies that the 
zeros of P w (z) are simple. Hence the images of these zeros under the above mapping are simple, 
except possibly when z = ±1/2. If P w (z) has a zero at z = ±1/2, then it follows that D w (p) may 
have multiple zeros at p = ±1. 

We now show that such multiple zeros at p = ±1 can only occur for a single value of w. The 
derivative of D w (p) at p = ±1 is given by 

D'Ml) = ±4(AV - 1) ± 2(A 2 - l)(p 2 - l)w. (3.6) 

If (A,//) (1, 1) then this derivative is zero when 

9 A V - 1 

" = V-ix^-i)- (3 - 7) 

Hence it is only at this single value of w that D w (p) can have multiple zeros at p = ±1. 

□ 

Under the transformation (|3.5p . the zeros of D w (p) are related to the zeros of the orthogonal 
polynomial P w (z), which are all real. The transformation from z to p then implies that D w {p) can 
only have roots on the unit circle or the real axis. Thus it makes sense to introduce the parametrisation 
p = e' lt for roots on the unit circle. A straightforward calculation leads to the results summarised in 
the next lemma. 
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Lemma 2. The zeros of the polynomial D w (p) are given by pk = e 



* , where 



tan wtk 




(3.8) 



( 



A 2 -n ( i? - 1 \ 



X 2 + l){^ + l) 



Additionally, for A = 1 or fj, = 1, D w (p) = atp = ±1. If A = fi = 1, p = ±1 is a triple zero of D w (p). 
In particular, for w sufficiently large, the numbers of zeros of D w (p), including their multiplicity, are 
given as follows: 

1. If A < 1 and [i < 1, equation \3. 8\) has 2w + 4 real solutions. D w (p) has 2w + 4 zeros on the 
unit circle. 

2. On the line segments A < 1 and [i = \, respectively A = 1 and \x < 1, equation hS. £)) has 2w + 2 
solutions. D w (p) has 2w + 4 zeros on the unit circle. 

3. If A = fi = 1, equation \3. 8\) has 2w solutions. D w (p) has 2w + 4 zeros on the unit circle. 

4- If A > 1 and fj, < 1 or A < 1 and /i > 1, equation §3. 8\) has 2w solutions. D w (p) has 2w zeros 
on the unit circle and 4 zeros on the real line. 

5. On the line segments A > 1 and \x = \, respectively A = 1 and fi > 1, equation 113. 8\) has 2w — 2 
solutions. D w {p) has 2w zeros on the unit circle and 4 zeros on the real line. 

6. If A > 1 and fj, > 1, equation \3. 8\) has 2w — 4 solutions. D w (p) has 2w — 4 zeros on the unit 
circle and 8 zeros on the real line. 

This corresponds to (2w + 4) — 4<t zeros p^ of D w {p) on the unit circle, where a S {0, 1,2}, depending 
on the values of X and [i. The other 4a zeros pt are located on the real line. If Pk is a real zero, then 
so is -p k , l/pk, and -1/pk- 

Proof. Equation (|3.8|) is obtained from D w (p) by substitution of p = e lt , followed by routine simpli- 
fication. The number of real solutions is most easily obtained by considering the graphs of the LHS 
and RHS of (j3.8H over the interval [0, 2ir). The function t&n wtk is monotonic and has 2w simple poles, 
whereas the behaviour of the RHS depends on the values of A and fj,. For example, if both A < 1 
and [i < 1, the RHS is monotonically decreasing and has 4 poles, leading to a total number of 2w + 4 
intersections of both graphs. The other cases can be obtained similarly. 



□ 



Remark. The real roots pk can be obtained correspondingly from pk = e Sk , where now 




(3.9) 
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4 Partition function identities 



In this section we derive explicit expressions for the partition function which are especially suited for 
an asymptotic analysis of the finite-size scaling behaviour. The key is the following Lemma. 

Lemma 3. 

Z„, w (a,b) = --^Res(/;pjfc) , (4.1) 

Pk 

where 

D w (p) V P J P 

and pk are the zeros of D w (p) . 

Proof. From equation f)2.8|> it follows that 

Z n , w {a,b) = J- I f(p)dp , (4.3) 



2vri _ 

with f(p) given in equation (|4.2|) , The integrand is a rational function in p so that the sum over all 
its residues on the Riemann sphere is equal to zero. From this it follows that 

Z n ^(a,6)=Res(/;0) = -^Res(/;p fc )-Res(/;cx)) , (4.4) 

Pk 

where the pk are the zeros of D w (p) . 
Due to the symmetry 

f(l/p) = -p 2 f(p) (4.5) 
the singularities p k come in pairs (p k , 1/pk) and the paired residues are equal, i.e. 

Res (f; Pk ) =Res(/;l/p fc ) . (4.6) 

It therefore follows that 

Z n , w (a,b) =Res(/;0) = -^Res (f; Pk ) . (4.7) 



2 

Pk 



□ 



It is convenient to first look at the special cases already considered above (see equation (|3.2p ). 
Proposition 4. We have 

2n+l w+i , , 

Z nw (l,l) = > sin 2 cos n , (4.8a) 

' v ' w + 2^ w + 2 w + 2' vy 

k=l 

r, / n 2 ™ 2 (2^ + l)vr _(2A; + l)7r , jol . 

u; + 1 2w + 2 2w + 2 

k=0 

2^+1 , , 

« 1| = iM^°" (4 - 8c) 

v ' k=0 

2 n ^ kir 

Z njW (2, 2) = — V cos™ — , (4.8d) 

fc=0 



S 



and 



i _ \2 \2w 

Z n , w (l + A 2 , 1 + A" 2 ) = i-f-A- (A + A" 1 )" (1 + (-1)") 



2A 2 1-A 2 " 



2 w COS" 



(^) 2 +sin 2 ^ 
for A / 1. 

Proof. Equation (|3.2p allows for an explicit calculation of the zeros of D w (p). An application of 
Lemma [3] gives the result. 

□ 

We note that the expression for Z n>w (l, 1) is a special case of results in [8]. For general values of 
A and n, we arrive at the following result. 

Proposition 5. For A/i ^ 1 and w sufficiently large, 

ZnMa ' ~ 4 ^ (A 2 - p 2 )(l - A 2 p 2 ) V^T J ^T^' (49) 



where p k are the zeros of D w (p) and 



A 2 1 /u 2 

+ 



i-\ 2 pI \ 2 -p\ i-fi> 2 pl v 2 -v\ 

Proof. Using Lemma El we write f(p) = r(p) / D w (p) with 



(4.10) 



r(p) = (1 - /)[(! - M V) + (/i 2 " pVH ( V i . (4.11) 



P J P 

According to LemmalU all roots of D w (p) are simple for w sufficiently large, and p = ±1 are removable 
singularities of f(p) due to the occurrence of the factor (1 — p 2 ) in r(p). Using that 

Res(/;p fc ) = r(p fc )/D( (p fc ) (4.12) 

for simple poles at pk, we find that 

Res(f;p k ) 2t«(A 2 -^)(m 2 -^)^ ,0 + 2^[(A 2 + m 2 -2A 2 m 2 ^)-(A 2 + M 2 -2^)^ u '] ' (4 ' 13) 
Eliminating p 2w using D w (p k ) = gives 

Res (f; Pk = — ^ — — . 4.14 

2 (A 2 -p|)(l - \ 2 p%) V Pfc / w + e k 

where e k is given by (|4.10p . □ 

Remark. The condition A/i 7^ 1 is necessary because when Afi = la root, p k , is located at p k = A, 
making the expression for Z ntW (a,b) invalid as it is written. However, in this case we already have 
given a much more explicit expression for Z n u ,(a,6) above. 



9 



5 Asymptotics 



If we consider the scaling behaviour of long walks in wide strips, the behaviour of our model should 
depend on the relation between average vertical displacement of an unrestricted polymer and the 
width of the strip, w. We therefore expect the occurrence of the scaling combination y/n/w in the 
asymptotics, and this is indeed what a mathematical analysis shows. 

As in the previous sections, it is convenient to first consider the special cases. 

Proposition 6. For n even, y/n/w fixed and n — ► oo, 

2 n 

Zn,wO-, 1) ~ ^3/5fao(Vn/ W ) ' ( 5 - la ) 

2 n 

Z ntW (l, 2) ~ — / 0>1 (^/ti;) , (5.1b) 
2 n 

Z ntW (2, 1) ~ —f lfi (rfi/ w ) , (5.1c) 

Z n , w (2,2) ~ ^f hl (^/ w ) , (5.1d) 



for A < 1 we have 



1 — A 2 2 n 

z n , w (i + A 2 , i + a- 2 ) ~ — ^ — ^ 2w ( A + A ~T + tfph^iyfrM > ( 5 - le ) 



while for A > 1 we have 



A 2 — 1 2 n 

Z n , w (1 + A 2 , 1 + A- 2 ) ~ (A + A" 1 )" + —f^rfi/w). ( 5 .lf) 

T/ie functions f\ jfi are given by 

/o l0 (x) = 2ttV f; fcV*^ = 27rV e -^ ^ (e"^) , (5.2a) 

O ■} V"" > - , 7T 2 (fc+l/2) 2 2 „ „ A 2 , / A 2 \ 

/o,i(x) = 2ttV (jfe + l/2) 2 e 2 x = 2ir 2 x 3 e — ~ tf' 2 ( e — ~ J , (5.2b) 

« / \ A _ * 2 (fc+l/2) 2 3 / _A£\ 

/i,o(x)=x 2^ e 2 =x-d 2 [e 2 j, (5.2c) 

= x f; e"^- 2 = x tf 3 (e-^) , (5.2d) 

/or A + I, 

(A 2 + l) n o A 2 2 (A 2 + l) o , A 2 , / A 2 \ 

A, A -i(x) = ^^2^V^fc 2 e-^ = ^^vrVe— r" tf' 3 (e~ *"J , (5.2e) 
#ere, 2 (?) = En=-oc 9 (n+1/2)2 and 3 (g) = E^=-oo ^ « re e// ^c ^-functions. 
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Proof. We shall only discuss the case A = fj, = 1, as the argument applies mutatis mutandum to the 
other cases. Upon introducing the variable x = n 1 / 2 /w, we can write 

fc=0 



For fixed k we find that 



lim cos n I — -pr I = e~^~ x , (5.4) 
n— >oo \Tl ' J 



and a similar contribution comes from the upper boundary of the summation for k 1 = w — k fixed. 
The contribution to the sum from other terms can be easily shown to be negligible, so that 

1/2 00 o o 

n ' , , v-^ A 2 „,2 

fc=— oo 

□ 

In Figure H] we have plotted the scaling functions at the points A,/j G {0, 1}. These show that the 
partition functions converge rapidly to the scaling functions even at quite modest widths (w ranges 
from 16 to 80 and the maximum length is lQw 2 ). The non-monotonic behaviour of foi(x) has also 
been observed in the three dimensional self-avoiding walk model at a comparable point in the phase 
diagram [3]. 

Remark. We note that even though we have only considered the scaling limit of x = \fnjw fixed, 
the above scaling forms permit asymptotic matching to the limits w — > oo with n fixed, and n — > oo 
with w fixed. For example, we obtain 



2 



J-2 n n -> oo 



—f x ^/w) ~ <^ (5.6) 



nl/2 ' I -^2 n w^oo 



which agrees with the known asymptotic behaviour of Z ntW (2, 2) for fixed w or n, respectively. 



In order to discuss the case of general a and b, let us consider equation (14.9)) more closely. Note 
that for large w the terms given by (|4.10p become negligible in comparison with w if Pk is a root 
on the unit circle, as is then uniformly bounded in k, w and n. Also note that in this case the 
dependence of (|4.9|) on [i in enters only through the roots pk- 

However, if pk is a real root it needs to be treated individually, and we therefore split the partition 
function into two sums, 

Z n!W {a, b) = Z(%(a, b) + Z(%(a, b) , (5.7) 

where 

^)fq^- 1 + A2 V ( 1 + p2 A n 1 (5 8) 



ii 



0.5 1 1.5 



scaling variable 



0.5 1 1.5 



scaling variable 



(/i,A) = (0,0) 



(/x,A) = (1,0) 




scaling variable 



scaling variable 



Ou,A) = (0,1) 



0M) = (1,1) 



Figure 4: T/ie scaling functions (solid lines) at (//, A) = (0, 0), (1,0), (0, 1) and (1, 1) respectively. Also 
plotted is numerical data (dashed lines) for widths ranging from 16 to 80; these show that the partition 
functions converge rapidly to the scaling forms. 



and 



4iM) 



l + A 2 



E 



2\2 



(i-p£) 



Ipfc I — J- 



2 \ n 



+Pk \ 



(5.9) 



Here, Zn}v(a, b) is a sum over zero, four, or eight terms, depending on the value of w, A, and /x. While 
it turns out that Znjv(d, b) admits a scaling form, Zn} w {a, b) does not. The following propositions give 
asymptotic estimates for Zn}u(a,b) and Z$ w (a,b) for large w. 

The analysis of Zn} u {a, b) uses an asymptotic estimate of the zeros for which p*. > 1, analogous to 
equations (7.7) and (7.10) in [2). If A > 1 and \i ^ A, we find 



p 2 ~ A 2 + (A 4 - 1) 



and if /J > 1 and A 7^ //, we find 



A 2 - /i 2 



A 



-2w 



(5.10) 



p 2^ M 2 + (M 4_ 1} 



2, ,2 



A Z /i 



/x 2 - A 2 



-2io 



(5.11) 
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while for A = fi > 1 we find 

p 2 ~ A 2 ± (A 4 - 1)A~ 



(5.12) 



I 1 i 


rr(l) 




rr(2) 
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A, 



Figure 5: Regions that have different scaling forms. The definitions via inequalities occur in equa- 
tion \5.14\ )- The region includes the dashed line segment A = 1,/z > 1 and the region 
includes the dashed line segment fj, = 1, A > 1 



In Figure [5] we introduce a finer division of the regions of the phase plane which will be needed 
below; essentially we split the regions T and B into two parts. 

From the considerations above we arrive at the following result. 



Proposition 7. For n even and n — ► oo, 



' (A 2 + - 1)V - 1) 
M 2 (m 2 - A 2 ) 2 



where 



I 1 W (» + ^) 



(A 2 + 1)( M 2 -1) 2 (^-1) 



-2tu 



A 2 -l 



>+a* T+-j^-(A+A 1 



A 2 1 

A 2 -l 
A 2 



A + A- 1 ) 



(A + A- 1 ) + 



(A 2 + l)(^ 2 -l) 2 (^-l) 

M 2(^2 _ A 2)2 



1\" 



A 2 -l 
A 2 



(A + A- 1 ) 



1\« 



Stb = {(^A)|/x = A>1} 

r( 1 ) = {(^,A)|/i>i>A} 

= {(^A) | A > 1 > fi} 



T( 2 ) = {( M ,A)| / u>A>l} 
£( 2 ) = {(/i, A) | A > /i > 1}. 



in i/te region 

in £/ie region 

on the line Stb, 

in the region B^ 2 \ and 

in the region B^ . 

(5.13) 

(5.14) 
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The asymptotic estimate for Zn} w {a,b) is done in two stages. First, we present an estimate for 
large w which is uniform in n. 



Proposition 8. 

41 M) 



A + A" 
4A 



1 \ nn 



W 



sin 2 if. 



x-x- 



cos n t k [l + Oiw- 1 )] (5.15) 



+ sin tk 



uniformly in n, where if. are the roots of h3. 8\) in [0, 7r). 



Proof. Apart from substituting = e ltk and simplifying, the main effort lies in obtaining an estimate 
for £fc. For A,/i ^ 1, we can estimate it directly from (|4.10|) 



. . A 2 + l ^ + 1 

M < — — - + 



If A = 1, we find that 



|A 2 -1| | M 2 -1| 
/i 2 + l 



\£k\ < 



|/i 2 - 1| 



(5.16) 

(5.17) 
□ 



and analogously for fi = 1. Finally, if both A = 1 and \i = 1 then = 0. 

For the scaling behaviour of Zn}u(a,b), we find identical scaling behaviour as in the special cases 
discussed above (except for A-dependent pre- factors). 

Proposition 9. For n even, y/n/w fixed and n — ► oo, 

( 2 n 



~ < 



n 



n 



jj^f{Vn/w) 



where 



x i? 3 2 
A 2 + l 



(A 2 - l) 2 
x i?2 I e 2 
A 2 + l 



2?r 2 x 3 e — — t?' 2 | e — ~ 



2vr 2 x 3 e — a - e - 



/or A + I, 
for\ = 1. 

/or A = 1 and fi = 1. 
/or A ^ 1 and = 1, 
/or A = 1 and jtt^ 1, and 
/or A 7^ 1 and // / 1. 



(5.18) 



(5.19) 



l(A 2 -l) 2 "" " " 3 
i/ere, #2(9) = a*"^ 1 / 2 ** 2 and ^3(9) = Xl^-oo ^ are e ^%dic d-functions. 

Proof. The derivation is done in complete analogy to the special cases discussed above. The only 
additional consideration is that now tk is only known asymptotically for w large. If either A 7^ 1 and 
/U7^1orA = )U = l, then for large w we find 



, 7r 

t k ~ K — 
W 



(5.20) 
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and otherwise (i.e. if A = 1 or fi = 1 but not A = /i = 1) for large w we find 

i fc ~(fc + l/2)-. (5.21) 

w 

□ 

We point out that these functions, which are effectively multiples of the scaling functions (see 
below), are elliptic t9- functions. Further they are independent of A and fj,, apart from that overall 
multiplicative factor, which is dependent only of A (as it depends on the half-plane limit). Finally we 
note that the result for A = \i = is equivalent to asymptotic results of Flajolet et al. [9] concerning 
the distribution of heights of binary trees since there is a correspondence between directed paths and 
binary trees. 



6 Discussion 



We now return to the issue of how the scaling forms calculated above interpolate between the half-plane 
and the infinite-slit. For the half-plane we have (recalling a = 1 + A 2 and b = 1 + fi 2 ) 



lim Z n , w (a,b) 



~ < 



2(1+A 2 ) 
(1-A 2 ) 2 



2™ 



2 n 

^172 



(A 2 -l) 
A 2 



A + A 



0<a<2(0<A<l) 
a = 2 (A = 1), 
a > 2 (A > 1), 



(6.1) 



as n — ► oo. Note that the scaling is independent of b, as it should be for w > n. For any finite w we 
expect that 

Z n ,w(a, b) ~ B w (a, b) fi w (a, b) n (6.2) 



where as w — > oo we have 



H w (a,b) 



-1> 



< a,b < 2, 

a > b and a > 2, 

a < b and 6 > 2. 



(6.3) 



If we denote the right-hand side of equation (|6.ip as S hp (n) for each value of a then a canonical 
general scaling Ansatz is 

Z n<w (a, b) ~ S^(n) / phasc . (6.4) 

The scaling function /phase depends on which phase or phase boundary of the infinite slit the values of 
a and b correspond (so as to match with the infinite slit phases). The exponent u±, on the other hand, 
should depend on the phase of the half-plane problem. 

We immediately note that there is something unusual here and that while v±_ is expected to be 
1/2 for < a < 2 it is expected to be for a > 2 ! In Martin et al. [3j it was proposed that a scaling 
theory like that in equation (16. 4|) (see equations (3.5), (3.6) and (3.7) of that paper) should only hold 
in the desorbed region of the infinite slit or on its boundaries. That is, we should expect a scaling 
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theory only when < a, b < 2 — this is precisely when v± = 1/2 > 0. This is what we have found for 
the directed walk case. 

Now, to match the scaling in equation (|6.ip (half-plane matching) one would require that 

/phase 0*0 ~ 1 as x -> . (6.5) 

This is a similar requirement to that of equation (3.6) of Martin et al. 0]. On the other hand the 
scaling of /phase as x — > oo needs to be considered in order to match the scaling of the infinite slit 
(equations (16. 2p and (16. 3|) ). 

Our results show that the partition function Z n ^ w (a, b) can be found as two parts and as 
given in equations (15. 13ft and (15.180 respectively. Importantly, for < a, b < 2 we have Z^ = and 
the form (|5.18p for Z^ is precisely that of (|6.4p with the half-plane matching behaviour (|6.5p obeyed 
by (I5.19p . The functions (|5.19p also obey the required infinite slit matching behaviour: that is, 



/phase (-c) < 



for a = b = 2 
x 3 e -^z 2 /8 for a < 2 and 6 = 2, 



xe k 2 x 2 /8 for a = 2 and b < 2, and 



(6.6) 



x 3 e -^x 2 /2 for a < 2 and 6 < 2 

as x — > oo. This large x behaviour allows for the recovery of (|6.2|) with (|6.3p . They also give us the 
first correction to the free energy as a function of w (see equations (7.3) and (7.6) of Brak et al. [2]). 
From equation (3.7) of Martin et al. [I] we find that the above results also agree with that generic 
prediction. 

For any value of (a, b) such that either a > 2 and/or b > 2, Z^ ^ and furthermore Z^ >> Z^ 
as n — > oo. As far as we can calculate, Z^ cannot be written in the scaling function form (|6.4|) . We 
note that the scaling form of Z^ acts as a correction to scaling to that of Z^ which allows for a 
matching of the total form Z^ + Z^ with the half-plane and infinite slit limits. Since the dominant 
part of the form of Z^ in (|5,13p is that of (|6,2p with the connective constant as in (|6.3p the infinite 
slit behaviour is clearly adhered to (although the coefficient B w {a, b) may go to zero in this limit). 

Consider equation (|5.13p . For b > 2 and a < 2 (region T«) Z<«0 ~ fi 2w ^ as w — > oo (n fixed 
and large) so for the half-plane one is left with Z^ with x — > 0. For a > 2 and b < 2 (region B^') 
the scaling form in (|5.13|) coincides with that of (|6.1|) for a > 2 so that half-plane limit is simple. For 
a = b > 2 (line Stb) this is also true. For a > 2 and b > 2 with a ^ b (regions and B^) there are 
two parts of the scaling form of Z^ with one part going to zero for w — > oo leaving the appropriate 
scaling corresponding to the half-plane. 

The above shows how complicated, though mathematically complete, the scaling picture can be 
for this problem. To summarise, we have calculated the scaling for large widths and large lengths 
of directed walks confined between two walls that interact with the walk. We explicitly demonstrate 
that the conjectured scaling theory [1] for polymers confined in such a manner holds exactly for this 
model. This theory holds when the polymer is in a desorbed state, or on the boundaries of this region 
in the parameter space, that is critically adsorbing. 



16 



Acknowledgments 

Financial support from NSERC of Canada, the Australian Research Council and the Centre of Excel- 
lence for Mathematics and Statistics of Complex Systems is gratefully acknowledged by the authors. 

References 

[1] DiMarzio E A and Rubin R J 1971 J. Chem. Phys. 55 4318-4336 

[2] Brak R, Owczarek A L, Rechnitzer A and Whittington S 2005 J. Phys. A 38 4309-4325 

[3] Janse van Rensburg E J, Orlandini E, Owczarek A L, Rechnitzer A and Whittington S G 2005 
J. Phys. A 38, L823-L828 

[4] Martin R, Orlandini E, Owczarek A L, Rechnitzer A and Whittington S G 2007 J. Phys. A 40, 
7509-7521 

[5] Janse van Rensburg E J, Orlandini E, and Whittington S 2006 J. Phys. A 39, 13869 

[6] Brak R, Essam J, Osborn J, Owczarek A L and Rechnitzer A 2006 J. Phys: Conf. Ser. 42 47-58 

[7] Andrews G E, Askey R and Roy R 1999 Special Functions in Encyclopedia of Mathematics and 
its Applications (Cambridge: Cambridge University Press) 

[8] Krattenthaler C, Guttmann A J and Viennot X G 2003 J. Stat. Phys. 110, 1069-1086 

[9] Flajolet P, Gao Z, Odlyzko A, and Richmond B 1993, In Combinatorics, Probability, and Com- 
puting, vol 2, 145-156. 



17 



